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Plants are sensitive to thermal and electrical effects; yet the coupling of both, known as thermo¬ 
electricity, and its quantitative measurement in vegetal systems never were reported. We recorded 
the thermoelectric response of bean sprouts under various thermal conditions and stress. The ob¬ 
tained experimental data unambiguously demonstrate that a temperature difference between the 
roots and the leaves of a bean sprout induces a thermoelectric voltage between these two points. 

Basing our analysis of the data on the force-flux formalism of linear response theory, we found that 
the strength of the vegetal equivalent to the thermoelectric coupling is one order of magnitude larger 
than that in the best thermoelectric materials. Experimental data also show the importance of the 
thermal stress variation rate in the plant’s electrophysiological response. Therefore, thermoelectric 
effects are sufficiently important to partake in the complex and intertwined processes of energy and 
matter transport within plants. 


I. INTRODUCTION 

Plants are complex biological systems, which feed, de¬ 
velop and function thanks to a variety of elementary 
and cooperative processes that ensure transport of en¬ 
ergy and matter in the form of mineral elements, carbo¬ 
hydrates, and hormones in watery solutions called sap. 
Two categories of sap are transported in conductive tis¬ 
sues: xylem sap transported in xylem, a continuum of 
dead cells forming a pipe that opposes little resistance to 
the transport of water and mineral elements from roots 
to leaves; and phloem sap, containing water and sugar, 
formed in the leaves by photosynthesis, and transported 
in phloem (sieve tube elements) from the leaves to the 
roots. Explanations for xylem sap transport includes 
the cohesion-tension theoryii^ and multiforce theoriesj^^— 
and the pressure flow hypothesis proposes a mechanism 
for phloem sap transporti^ Notwithstanding their rela¬ 
tive merits, it is of interest to note that as saps contain 
ions, a net voltage may be recorded within plants. Fur¬ 
thermore, as a living plant is naturally submitted to a 
temperature gradient between its roots and its leaves, 
heat may also flow from its hotter to its colder regions. 
These two facts lead to the questions of the existence and 
strength of a coupling between electrical current and heat 
flux in plants, and whether this coupling plays a signifi¬ 
cant role in the plants’ lives. Answers to these questions 
are of importance to gain further insight into plant home¬ 
ostasis, closely related to the Le Chatelier-Braun princi¬ 
ple in chemical thermodynamics, and to develop further 
the field of plant electrophysiology by analysing the im¬ 


pact of thermoelectric coupling on photosynthesis, plant 
respiration, and any other processes involving bioelectro¬ 
chemical phenomena»i 

Biophysical processes, which include energy and mat¬ 
ter transport in dedicated structures, are essentially 
nonequilibrium in nature. Their description involves 
the notions of thermodynamic forces and their conjugate 
fluxes which, in a linear response description, are pro¬ 
portional to each other. In a general manner, transport 
phenomena are nonequilibrium irreversible processes;^ 
fluxes within a system result from applied external con¬ 
straints, and they are accompanied by energy dissipa¬ 
tion and entropy production. As the forces applied to 
a system derive from potentials, the description of the 
system’s nonequilibrium properties relies on the defini¬ 
tion of these potentials and their degree of couplingi^S 
Thermoelectricity is sometimes mistakenly viewed as a 
phenomenon pertaining to solid-state systems /^i^^ but 
this effect has also been reported in liquidsi^ and gels»i^ 
Now, considering more specifically electrical phenomena 
in plants, one should also account for their dependence 
on temperature and the local variations of this latter. As 
a matter of fact, one learns from the equilibrium ther¬ 
modynamics of charged fluids that heat and electricity 
are fundamentally related through the Gibbs-Duhem re¬ 
lation; so, since plants, as any other system, are subjected 
to the laws of thermodynamics, one may anticipate that 
a proper understanding of the coupled transport of heat 
and electricity in plants will pave the way towards the 
potentially rich new field of thermoelectricity-based plant 
energy conversion, and harvesting metabolic energy from 
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plants 

Thermoelectricity as part of non-isothermal processes 
in living organisms, has already been mentioned; but ex¬ 
cept for the case of insects like hornets and beesj^^— 
for which genuine thermoeletric effects have been anal¬ 
ysed, the other studies simply reported the record of 
the temperature rise of Colocasia odora leaves using a 
thermocouple?^ or analysed pyroelectricity, mistakenly 
taken as a thermoelectricity, amongst the nonequilibrium 
processes at the origin of life.— We propose here to con¬ 
sider bean sprouts as an illustrative case of a living plant 
system whose physiological properties also entail ther¬ 
moelectric effects, with a particular focus on the thermo¬ 
voltage response of the plant. While a voltage difference 
along the stem of a plant has already been measured 
and analysed)^ its interpretation does not account for a 
likely the thermoelectric origin. Nevertheless, it is clear 
that the source of any voltage difference is due to concen¬ 
tration gradients of charges inside the plant, regardless 
of the cause of the gradients. The underlying processes 
may thus originate in various electrochemical processes, 
including acid-base and redox effects— In the present 
work we focus on the voltage response due to a temper¬ 
ature difference. The system is presented in Fig. 1. The 
stem and leaves are kept at constant temperature while 
the temperature of the roots can be modified; the ther¬ 
moelectric response is recorded on both the plant (a) and 
a control wire (b). 



II. THEORETICAL CONSIDERATIONS 

Thermoelectricity differs from electrochemistry as the 
latter considers isothermal systems only. Thermoelec¬ 
tric effects, on the contrary, manifest themselves in non- 
isothermal conditions and their thermodynamic descrip¬ 
tion is best done with the coupled variables (^, T), where 
/i is the electrochemical potential and T the temperature 
of the considered system. In a force-flux approach, as a 
temperature bias is applied to a conductor, both heat and 
charges are transported in coupled flows since each elec¬ 
tron carries an electric charge, e, and energy; the ther¬ 
moelectric coupling parameter or Seebeck coefhcient a is 
thus defined as: 


= ( 1 ) 

e 

where ^ denotes the spatial gradient. Equation m may 
be viewed as a generalization of the Gibbs-Duhem rela¬ 
tion to the dynamical response of the coupled intensive 
variables by defining the thermodynamic forces 

acting on the system as the gradients of the thermody¬ 
namic potentials. Each force is conjugated to a flux which 
is proportional to the time derivative of the correspond¬ 
ing extensive parameter, and it follows that forces and 
fluxes are linearly combined as proposed by Onsager.— 
Note that the coefhcient a, which in Eq. © is dehned lo¬ 
cally, is also often expressed as a ^ —eAV/AT, where A 
denotes the global difference between two distant points 
of the system. 

These general considerations apply for any system con¬ 
taining free charge carriers. As a consequence, for a given 
temperature difference, the measurement of a thermo¬ 
voltage directly gives an estimation of the average ratio 
of charge concentrations in separate parts of the system. 
Therefore, any system containing free mobile carriers 
may be exhibit a thermoelectric signature. One should 
note that depending on the specihc properties of the con¬ 
sidered system, this signature may intertwine or not with 
those of other temperature-dependent processes; hence 
the need for specihc system-dependent approaches to ex¬ 
tract the thermoelectric signal. In the case of liquids and 
gels, the studies usually focus on the thermoelectrophore¬ 
sis parameters— For solid state-systems, the con¬ 
duction is ensured by mobile electrons and holes, leading 
to a possible steady-state electrical current. For ionic 
carriers, in liquids and gels, such a steady-state current 
may only take place if a redox process occurs at the ends 
of the system. 


FIG. 1. General view of the experimental setup. The voltage 
records are respectively on the plant (a), and the control wire 
(b). 


III. EXPERIMENT 

The experiment aims to identify and characterize the 
thermoelectric response of a living plant (a bean sprout 
in the present work). As the response of the plant nec¬ 
essarily entails physiologic effects, the measured data is 
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compared to the purely thermoelectric response of a con¬ 
trol wire. As depicted on Fig. [H the temperature of 
the plant’s roots is imposed by a thermostatic bath; two 
thermocouples are used to measure the temperatures of 
the roots and the leaves, and two electrodes are placed 
respectively in contact with the roots and the leaves 
of the liviirg plant. The roots are bathed in an aque¬ 
ous dilute KCl solution. The electrodes connected to 
the leaves are made of Ag/AgCl wires. All the tem¬ 
peratures and voltages are recorded using a Keithley 
K2700 scanning multimeter and a K182 nanovoltmeter. 
The thermoelectric control wire (CW) is made of two 
chromium-nickel-steel junctions, with a Seebeck coeffi¬ 
cient ajunction = 4.6 /K; it is placed in the same config¬ 

uration as that of the plant. The recorded voltage at the 
ends of the control wire describes a pure thermoelectric 
response, exempt from any influence of the physiologic 
process. 



beck coefficient. Conversely, the plant exhibits a com¬ 
plex response, which includes the effects of its physio¬ 
logical processes thus giving rise to an active Seebeck 
response. The curves represented on the lower panel are 
displayed together on the same scale to qualitatively show 
how the responses (voltages and Tpiant) follow the 
applied constraint (temperature difference AT). The fig¬ 
ures on the y-axis correspond to AT in °C. The values 
taken by the voltages Ccw and l^iant may be determined 
from the curves using a scaling factor that corresponds to 
the relevant Seebeck coefficient. The control wire shows 
a faithful homothetic response to the temperature dif¬ 
ference, with a linear factor of scaling of 4.6 yV/K as 
expected for the Seebeck coefficient of the correspond¬ 
ing chromium-nickel-steel thermocouple. The response 
of the plant gives an average Seebeck coefficient of 2.5 
mV/K. Such a large value in a solid-state material would 
correspond to a system with low electrical conductivity. 
Here, this results from ion movements through xylem and 
phloem saps circulating in opposite directions. 

To analyze precisely these data, it is convenient to split 
each panel of Fig. [5] into three parts, which correspond 
to three distinct time domains. The first and third ones, 
Figsl3]and[5l respectively, describes the plant’s response 
when its roots are heated; while the second, FigUl reports 
the response over the cooling sequence. 



FIG. 2. Experimental data including the voltage response 
of the control wire I4w (upper left panel) and of the plant 
Vpiant (upper right panel), both submitted to a time-varying 
temperature difference AT, reported on the lower panel. The 
magnitude of measured voltage I4w has been multiplied by 
a factor 2000 in order to report it on the same scale as that 
of Vpiant- The red-dotted lines indicate the average slope of 
the collected data. In the lower panel, the scale is that of the 
temperature difference; the voltages are also depicted in this 
lower panel to see how close their time-dependent behaviour 
is to that of the measured temperature difference. 

The voltage response of the living plant and that of 
the control wire to the temperature difference AT are 
both shown on Fig. [21 where the time evolution of AT 
is also displayed, illustrating the dynamics of the plant 
and that of the control wire. These data were recorded 
over a period of four hours for various temperature dif¬ 
ferences and variation rates. The two upper panels show 
separately the parametric plots characterising the ther¬ 
moelectric response of the control wire and that of the 
plant. Notice that the response of the control wire is 
strictly proportional to the temperature difference, as ex¬ 
pected for a thermoelectric material, with a passive See- 
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FIG. 3. First heating sequence. The plant’s response is 
strongly influenced by the temperature variation rate dT /dt. 


Depending on the precise thermal condition supported 
by the plant, we can identify different thermoelectric re¬ 
sponses. Consider the first warming sequence of Fig. |31 
One may identify five events (labelled from a to e). Dur¬ 
ing the first period (Fig. |3^), the temperature difference 
is roughly constant. The plant presents a large fluctu¬ 
ating thermoelectric response while the control wire is 
characterized by a constant Seebeck coefficient. During 
the second period, shown on Fig. [Slo, the plant’s reac¬ 
tion is essentially sensitive to the temperature variation 
rate dT/dt. This sensitivity has already been observed 
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FIG. 4. Cooling sequence. The response is roughly linear 
but visible discrepancies are due to the plant’s physiological 
processes. 



FIG. 5. Second heating sequence. 


in studies devoted to plant sensingj ^^i^^ The processes 
taking place within the zones a and b can be also ob¬ 
served in the c and d zones; and when the rate dT/dt 
is reduced, we observe a linear response for the plant, 
as shown in the zone e. It thus appears clearly that 
the passive response, probably mostly due to the xylem 
properties, can be strongly modified by active physio¬ 
logical processes. This latter appears to be particularly 
sensitive to the time derivative of the stress, but not so 
much on the intensity of the stress itself; this observation 
is perfectly consistent with previous findings 

We now turn to the cooling sequence. Fig. S] The 
cooling rate is approximately constant, except for the 
two events when it is clearly larger than the average as 
indicated with the labels a and b, on the upper right 
panel. These events relate to time delays, also apparent 
in the measured data for the control wire on the upper left 
panel: these are signatures of the change of rate dT/dt. 


Note that although these may serve as a probe, time 
delays must be considered with care, since in the case 
of very large time response of a system, they may lead 
to erroneous interpretations^^ The shape of the plant 
response for events taking place in the zones a and b is 
very similar to that of Figs. [3 }d and |3Jl, which conhrms 
the acute sensitivity of the plant to the rate dT/dt. This 
physiological response is the fastest response to a stress 
that a plant can use, as it has no means to avoid it. 

The last sequence depicted on Fig. [5] concerns the 
plant’s response to warming at a reduced rate. It is essen¬ 
tially linear, thus confirming the presence of a threshold 
rate process for the plant’s response to a given stress. It 
also shows unambiguously that thermoelectricity natu¬ 
rally partakes in the complex processes that contribute 
to a plant’s life. 


IV. DISCUSSION AND CONCLUDING 
REMARKS 


All the effects discussed in the present work are related 
to the internal modihcations of the ionic concentrations 
inside the plant, including the regulation of cellular ion 
transporters allowing ion release into the xylem sap to¬ 
ward the shoots. It is clear that while the response is 
mainly driven by a classical thermoelectric process with 
a giant Seebeck coefficient, the physiology of the plant 
interferes in a complex way with it. A complete nu¬ 
merical model, including series and parallel assemblies 
of living (phloem sieve tube elements, root and endoder- 
mis cells) and dead cells (xylem tracheids and vessels), 
will be considered for further investigation of the inter¬ 
play between physiologic and thermoelectric response to 
a thermal stress. Further, since there is no limitation 
in the number of thermodynamic potentials which may 
be coupled)^ we may account for pressure and derive a 
complete response of a system to the direct and the cross¬ 
effects of the three thermodynamic potentials, {P,^,T). 
Then the lowering of a local pressure gradient, by the 
coupled effect of the temperature and the electrochemi¬ 
cal potential becomes possible. This latter process may 
contribute to the motion of saps in plant, and especially 
in trees. 

The reported measurements are fully scalable from the 
unique cell size to a complete plant or tree. In addi¬ 
tion to the stationary response observed here, the voltage 
fluctuations are also of great interest since they permit 
the study of the threshold of the physiological response 
through its noisy signals. As observed in other systems, 
the noise response is a very sensitive probe of the emer¬ 
gence of a macroscopic response Taking each cell as 
a fluctuator, the convolution of the individual responses 
may lead to different signatures depending on the con¬ 
sidered scale of the sample. 
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Appendix A: Plant material culite particles and then transferred in dilute KCl solu¬ 

tion (120 mM) for experiments. 


Seeds of bush bean (Phaseolus vulgaris L. cv. con¬ 
tender) were sown into pots containing vermiculite as 
soil. The bean seedlings were grown in a growth cham¬ 
ber at 25 ± 2°C with a cycle of 12 hours of light (40 /xmol 
photon m~^- s“^), 12 hours of darkness. The plants were 
watered at the bottom of the pots every three days. The 
21 days old plants were removed from the vermiculite. 
Roots were carefully rinsed with water to remove vermi- 


Appendix B: Measurement electrodes 

Two electrodes were placed respectively at the surface 
of the root and the stem of the living plant. The elec¬ 
trodes connected to the stems are made of Ag/AgCl wires 
moistened with 200 mM KCl and wrapped in cellophane 
to provide appropriate contact with the plant surface. 
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